Study on Satellite-to-ground Laser Communications with Existence of Atmosphere

The laser beam used to establish a communication channel between satellite and ground segments has a small divergence angle and a tiny spot on the Earth’s surface, which may lead to the fail of the system. So it’s important to study the deflection of laser beam by the Earth’s atmosphere and find a way to calibrate this error. Both theoretical analysis and real data processing method are used to obtain the mathematical model for divergence angle of laser communication beam and its correction function. Then the model has been applied to the data, which was used to describe the atmosphere state by traditional ground segments to obtain the critical elevation angle. According to the results of calculation, our conclusion will be that the correction should be done when the critical elevation happens.      

基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

导航卫星中断频次的分析方法

面向导航卫星中断频次的定量分析需求,该文分析导航卫星中断产生的主要原因,给出中断频次分析的流程,并针对分析过程中的三个关键问题,研究提出具体实施方法,包括通过相关性分析快速定位底层中断事件,通过中断树建立指标分析模型,并融合在轨数据、地面试验数据快速预估得到底层功能异常率等。最后通过示例进一步说明中断频次分析过程。该文方法已应用于北斗导航卫星工程。     

抗辐照通用扩展SiP芯片设计

研究微小卫星综合电子系统的SiP技术实现方法。首先介绍微小卫星综合电子系统结构的组成和采用SiP技术的必要性，然后对综合电子系统进行功能模块划分，并对其通用扩展模块进行详细的SiP设计，包括抗辐照器件选型、原型验证、SiP原理图、基板管壳一体化设计、建模仿真、制造加工、实装测试验证等，通过SiP技术实现了一种星载综合电子系统中通用扩展SiP芯片产品，经过实际验证测试，在保证模块功能和性能的前提下，整体模块重量从230 g减轻到48 g，体积由180 mm?130 mm?17 mm减小到46 mm?46 mm?8 mm，很好地满足了星载功能模块小型化、轻量化设计需求。     

星载激光通信终端在轨机动对温度影响

对激光通信终端在轨瞬态温度变化开展了仿真，以期研究在轨机动的影响、热分析和热试验时机动模式的简化模拟。通过合理地分析与简化，建立终端的轨道热分析模型，准确模拟在轨机动，根据典型机动模式、外热流和涂层退化等因素设计了计算工况，得出了终端在轨运行过程中的温度场随时间和姿态变化规律。结果表明：不同机动模式下的温度变化存在差异，最大可达23.0℃，采用固定姿态或一维转动的简化热分析或热试验不能准确模拟实际飞行温度，甚至不能部分替代二维转动热分析或热试验；光学天线和反射镜的温度控制是制约终端工作的瓶颈，为终端设计合适的避光机动策略，可大幅度提高温度稳定度和均匀度。研究结果可以为光机电设备在轨机动策略的设计和改进提供参考。     

Fully distributed time-varying formation tracking control for multiple quadrotor vehicles via finite-time convergent extended state observer

This paper investigates a time-varying anti-disturbance formation problem for a group of quadrotor aircrafts with time-varying uncertainties and a directed interaction topology. A novel Finite-Time Convergent Extended State Observer (FTCESO) based fully-distributed formation control scheme is proposed to enhance the disturbance rejection and the formation tracking performances for networked quadrotors. By adopting the hierarchical control strategy, the multi-quadrotor system is separated into two subsystems: the outer-loop cooperative subsystem and the inner-loop attitude subsystem. In the outer-loop subsystem, with the estimation of disturbing forces and uncertain dynamics from FTCESOs, an adaptive consensus theory based cooperative controller is exploited to ensure the multiple quadrotors form and maintain a time-varying pattern relying only on the positions of the neighboring aircrafts. In the inner-loop subsystem, the desired attitude generated by the cooperative control law is stably tracked under a FTCESO-based attitude controller in a finite time. Based on a detailed algorithm to specify the cooperative control protocol, the feasibility condition to achieve the time-varying anti-disturbance formation tracking is derived and the rigorous analysis of the whole closed-loop multi-quadrotor system is given. Some numerical examples are conducted to intuitively demonstrate the effectiveness and the improvements of the proposed control framework.     

An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

Evaluation of a regional real-time precise positioning system based on GPS/BeiDou observations in Australia

The performance of real-time (RT) precise positioning can be improved by utilizing observations from multiple Global Navigation Satellite Systems (GNSS) instead of one particular system. Since the end of 2012, BeiDou, independently established by China, began to provide operational services for users in the Asia-Pacific regions. In this study, a regional RT precise positioning system is developed to evaluate the performance of GPS/BeiDou observations in Australia in providing high precision positioning services for users. Fixing three hourly updated satellite orbits, RT correction messages are generated and broadcasted by processing RT observation/navigation data streams from the national network of GNSS Continuously Operating Reference Stations in Australia (AUSCORS) at the server side. At the user side, RT PPP is realized by processing RT data streams and the RT correction messages received. RT clock offsets, for which the accuracy reached 0.07 and 0.28?ns for GPS and BeiDou, respectively, can be determined. Based on these corrections, an accuracy of 12.2, 30.0 and 45.6?cm in the North, East and Up directions was achieved for the BeiDou-only solution after 30 min while the GPS-only solution reached 5.1, 15.3 and 15.5?cm for the same components at the same time. A further improvement of 43.7, 36.9 and 45.0 percent in the three directions, respectively, was achieved for the combined GPS/BeiDou solution. After the initialization process, the North, East and Up positioning accuracies were 5.2, 8.1 and 17.8?cm, respectively, for the BeiDou-only solution, while 1.5, 3.0, and 4.7?cm for the GPS-only solution. However, we only noticed a 20.9% improvement in the East direction was obtained for the GPS/BeiDou solution, while no improvements in the other directions were detected. It is expected that such improvements may become bigger with the increasing accuracy of the BeiDou-only solution.